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ABSTRACT 

We investigate the influence of scattering and geometry on the attenuation curve in disc 
galaxies. We investigate both qualitatively and quantitatively which errors are made by 
either neglecting or approximating scattering, and which uncertainties are introduced 
due to a simplification of the star-dust geometry. We find that the magnitude of these 
errors depends on the inclination of the galaxy, and in particular that for face-on 
galaxies, the errors due to an improper treatment of scattering dominate those due to 
an imprecise star-dust geometry. Therefore we argue that in all methods which aim 
at determining the opacity of disc galaxies, scattering should be taken into account in 
a proper way. 
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1 INTRODUCTION 

During the last decade it has become clear that also in ex- 
tragalactic astronomy, dust obscuration plays an important 
role, and that the construction of accurate galaxy models 
■ requires a proper radiative transfer modelling. However, it 
' is often very time-consuming, and therefore practically un- 
, feasible, to solve radiative transfer problems exactly. For ex- 
' ample, assume we want to know the spatial distribution of 
stars and dust for a certain galaxy, where the only known 
, data are observed surface brightnesses. The only means to 
' solve this problem is to assume a certain parametrized model 
for the galaxy, containing a stellar and a dust component, 
and to solve the radiative transfer equation (RTE) for each 
set of parameters, and then select the model that fits the 
observations best. Such a modelling procedure demands an 
enormous number of RTE solutions, because the parameter 
space has to be sufficiently large in order to produce real- 
istic models, in particular for disc galaxies (Misiriotis et al. 
2000). Therefore it would be very useful if we could simplify 
the RTE in any way. 

A first way to simplify the RTE is rather obvious: 
neglect the scattering term, which is responsible for the 
integro-differential character of the RTE (Christensen 1990, 
Jansen et al. 1994, Ohta & Kodaira 1995). However, al- 
ready more than a decade ago, Bruzual et al. (1988) clearly 
demonstrate the importantance of including scattering in 
radiative transfer calculations. Several other authors have 
since then confirmed these warnings, in particular Witt et 
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al. (1992). Another option to simplify the RTE is an approx- 
imation of the scattering term, by either forward or approx- 
imate isotropic scattering, or a combination of both (Natta 
& Panagia 1984, Guiderdoni & Rocca-Volmerange 1987, 
Calzetti et al. 1994). Di Bartolomeo et al. (1995) compare 
the attenuation curves calculated with these various ap- 
proximate RTE solution methods, and find significant dif- 
ferences. Although these differences can partly be ascribed 
to different assumptions concerning the optical properties 
of the dust, it is very probable that also the approximate 
solution methods are responsible for a significant error. 

It is obviously important to know the errors introduced 
on the attenuation curve by the various approximate so- 
lution techniques, but these have never been clearly deter- 
mined. Either the approximate solutions are simply adopted 
without concern about the errors, or the errors are estimated 
for very simple galaxy models, such as homogeneous slabs 
or sandwich models (Bruzual et al. 1988, Di Bartolomeo 
et al. 1995). To make things worse, Disney et al. (1989) 
showed that the geometry of the stellar and dust distribution 
strongly determines the observed radiation field (without in- 
cluding the effects of scattering however). Lack of knowledge 
about the geometry can thus also introduce significant errors 
on the attenuation of galaxies. 

In a previous paper (Baes & Dejonghe 2001, paper I) 
we described four different methods to solve the RTE in a 
plane-parallel geometry. These methods accomodate an ar- 
bitrary vertical distribution of stars and dust. We will now 
use these methods will allow us to address both issues of 
the RTE mentioned above. On the one hand, we can adopt 
them for a disc galaxy model with a realistic vertical struc- 
ture, and compare the results with those obtained by us- 
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ing the various approximate solutions. This wiU allow us to 
quantify the errors introduced by the different approxima- 
tions, and the importance of properly including scattering 
effects. On the other hand, we can apply these methods to 
a family of realistic galaxy models that can accomodate a 
wide range in distribution of stars and dust. This will allow 
us to investigate the influence of the geometry of the stellar 
and dust components on the attenuation curve, without any 
simplifying assumptions on the RTE such as the neglect of 
scattering. Last but not least, wc can solve our RTE prob- 
lems with four completely different methods: consistency is 
then a guarantee for accuracy. 

In Section 2 wc describe the radiative transfer mecha- 
nism and the ways to obtain the solution, and in Section 3, 
we present our set of disc galaxy models. In Section 4 and 5 
we discuss respectively the influence of scattering and ge- 
ometry on the attenuation curve as described above. In Sec- 
tion 6 we discuss the results. 



2 RADIATIVE TRANSFER MODELS 
2.1 The RTE in plane-parallel geometry 

In plane-parallel geometry, the RTE can be written as 
dl 

+ \ujk{z) J'^ I(z, m') V, (1) 

where I{z,fi) is the specific intensity of the radiation at a 
height z above the plane of the galaxy, and in a direction 
which makes an angle arccos /x with the face-on direction 
^ = 1. The known quantities in this equation are the dust 
opacity the stellar cmissivity 77,(2), the dust albedo 

u> and the angular redistribution function (hereafter ARE) 
*(/i, /u'). The RTE can be brought in another form by in- 
troducing the optical depth r instead of z, 

t,{z')dz', (2) 

which yields 

M 1^ (r, m) = I{r, n) - S, (r) - | ^ /(r, mOV, (3) 

with S'*(r) = rit{r) / k{t) is the stellar source function. This 
equation is to be solved for ^ r ^ To, the total optical 
(face-on) depth of the galaxy, 

/oo 
K{z)dz. (4) 
-00 

Although the RTE can be solved for any point in the galaxy, 
we will focus on the attenuation A{jj), i.e. the fraction of 
the intensity that is attenuated by the dust detected by an 
observer at r = 0, into a certain direction /i. Because the 
RTE is a linear equation, this fraction will be independent 
of the total amount of stellar emission. We can thus choose 
the normalization of the stellar emissivity (or the source 
function). We take 

/oo /-TO 
r)*{z)Az= / 5,(r)dT = l, (5) 
-00 Jo 



which means that the intensity that leaves the galaxy in the 
absence of dust in the face-on direction equals 1. In another 
direction ;u, the dust-free intensity is then simply l//u and 
the attenuation (in magnitudes) is 

^(M) = -2.51og[M/(0,M)]. (6) 



2.2 Solution of the RTE 

Because of the complexity of the RTE, a lot of authors have 
tried their ingenuity to find sophisticated methods to solve 
it. In paper I we presented four different methods to solve the 
RTE in plane-parallel geometry, which can handle absorp- 
tion and multiple scattering, arbitrary vertical distributions 
of stars and dust and arbitrary phase functions. 

(i) The spherical harmonics method consists of expanding 
all the angle-dependent terms in the RTE into a series of 
spherical harmonics, such that the RTE is replaced by a set 
of ordinary differential equations. The method turns out to 

be very efficient. 

(ii) In the discretization method, originating from the 
theory of stellar atmospheres, integrals are replaced by sums 
and differentials by finite differences, resulting in a set of 
vector equations which can be solved iteratively. 

(iii) In the iteration method, the intensity is expanded in 
a series of partial intensities. Each of the partial intensities 
obeys its own RTE, which can be solved iteratively. It can 
easily be extended to more complex geometries, but com- 
pared to the spherical harmonics method it is a very costly 
algorithm. 

(iv) The Monte Carlo method is probably the most widely 
adopted method to solve radiative transfer problems, and 
the basic idea of the method consists of following the trajec- 
tory of a large amount of individual photons through the 
galaxy, whereby the fate of a photon on its path is de- 
termined by random events. For one-dimensional radiative 
transfer its efficiency is comparable to that of the iteration 
method. 

For more details about these methods and references to the 
literature we refer to paper I. In this paper we will study the 
radiative transfer through a plane-parallel slab, as a func- 
tion of the different physical and geometrical parameters. 
Having these four different methods to solve the RTE at our 
disposal proves to be very useful, not only to check the ac- 
curacy of the numerical results, as already mentioned, but 
also to understand the physical background of observed phe- 
nomena. 



3 THE DISC GALAXY MODELS 

In order to construct disc galaxy models we have to charac- 
terize the functions that appear in the RTE (1), i.e. specify 
the physical properties and spatial distribution of stars and 
dust. We will adopt a set of realistic galaxy models with a 
rnimbcr of parameters. Varying these parameters will then 
enable us to investigate the influence of scattering and geom- 
etry on the attenuation curve. We will also adopt a template 
model to compare other models with. 
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Table 1. The adopted data sets for the optical properties of the 
dust grains. Tabulated are the optical depth t relative to the V- 
band value, the scattering albedo ui and the asymmetry parameter 
9- 



band 


A (fim) 


r 




g 


UVl 


125 


3.44 


0.42 


0.60 


UV2 


150 


2.75 


0.42 


0.59 


UV3 


175 


2.44 


0.48 


0.58 


UV4 


200 


2.87 


0.48 


0.54 


UV5 


225 


3.01 


0.51 


0.46 


UV6 


250 


2.38 


0.57 


0.46 


UV7 


300 


1.92 


0.57 


0.47 


U 


360 


1.60 


0.57 


0.49 


B 


440 


1.32 


0.57 


0.48 


V 


550 


1.00 


0.57 


0.44 


R 


700 


0.73 


0.54 


0.37 


I 


850 


0.47 


0.51 


0.31 



3.1 The optical properties of the dust 

Although it has been shown that the physical properties of 
dust in different environments can vary greatly (Witt et al. 
1984, Mathis & Cardelli 1992), we will assume, for sake of 
simplicity, one single kind of dust grains. This means that 
the spatial and wavelength dependencies of all quantities 
appearing in (1) are separable. In particular, the dust albedo 
and the ARF are then independent of position. 

To describe general anisotropic (conservative) scat- 
tering we adopt Henyey-Greenstein scattering (Henyey & 
Greenstein 1941). The ARF corresponding to this kind of 
scattering is a one-parameter function parametrized by the 
asymmetry parameter g, which is the average of the cosine 
of the scattering angle. A closed expression and a plot of 
the Henyey-Greenstein ARF can be found in Appendix A 
of Paper I. In Section ^ other kinds of scatterings such 
as forward and isotropic scattering, will be compared with 
Henyey-Greenstein scattering. 

There are two ways to determine the wavelength de- 
pendence of the albedo u, the asymmetry parameter g and 
the total optical depth tq. On the one hand, values can be 
derived theoretically, by assuming a certain dust grain com- 
position and calculating the optical properties of the dust 
(e.g. Draine & Lee 1984). On the other hand, optical prop- 
erties can be derived empirically, usually based on a variety 
of observations of scattered light in the Galaxy (e.g. Bruzual 
et al. 1988). Our data set consists of the optical properties 
theoretically derived by Maccioni & Perinotto (1994), dis- 
played in Di Bartolomeo et al. (1995). The adopted values 
are tabulated in table |l] for the central wavelengths of vari- 
ous bands. 

3.2 The stellar distribution 

The vertical distribution of stars in disc galaxies is still a 
matter of debate. The most straightforward way to derive 
such a vertical distribution is to study the surface bright- 
ness of edge-on galaxies at different heights above the plane. 
There is a general consensus that at great heights the light 
distribution decreases exponentially. Close to the plane of 
the galaxy however, dust attenuation makes the observation 
of the vertical distribution difficult. As a consequence, differ- 



1 - 




z/h 

Figure 1. The stellar emissivity rif{z) of the family of sech' 
models (^. It is shown for different values of the parameter q, 
ranging between the extreme values g = and q = 2, which are 
indicated and plotted in bold. It is clear that they all have the 
same exponential asymptotic behaviour, whereas the sharpness 
of the peak at 2 = is determined by q. 



ent models have been proposed. The most popular models 
are an isothermal sheet distribution (van der Kruit & Searle 
1981) and an exponential profile (Wainscoat et al. 1989). 
In order to allow a wide range in vertical density distribu- 
tions, we adopt a two-parameter family of stellar emissivity 
profiles which were adapted from van der Kruit (1988), 

where the function B[x,y) in the represents the Beta func- 
tion (Abramowitz & Stegun 1972). The first factor is deter- 
mined such that the normalization condition (^) is satisfied. 
At great heights above the plane of the galaxy, the light 
distribution (H) indeed decreases exponentially, 

ri4z) oc e"'''!/'' for \z\ > h. (8) 

The first parameter in ^ is the exponential scaleheight h. 
It is generally known that the velocity dispersion of different 
stellar types in the Galaxy is correlated with the age of the 
stars (e.g. Dehnen & Binney 1998). Due to this dynamical 
heating of the disc, the vertical scaleheight of the stars in- 
creases as a function of age (Wainscoat et al. 1992). Because 
the younger stellar populations are bluer than older ones, it 
can be expected that the scaleheight of stars increases with 
wavelength. Xilouris et al. (1999) accurately modelled sev- 
eral edge-on disc galaxies and found no significant trend of 
the stellar scaleheight with wavelength (from the blue to 
the near-infrared K band). Therefore we adopt a constant 
stellar scaleheight. 

The second parameter in (^, g, corresponds to the 2/n 
used by van der Kruit (1988), and takes values between 
and 2. It determines the shape of the stellar density near the 
plane of the galaxy, and hence is called the shape parameter. 
More precisely, the peak of the emissivity at a = becomes 
sharper as q decreases. A number of emissivity profiles for 
various values of q is illustrated in Figure |^. For the special 
cases corresponding to the values g = 0, 1 and 2, we find the 
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exponential, the sech- and the isothermal disc respectively, 



els.^ It is fairly straightforward to check that 5, (r) will have 
the form 



g = 
g=l 
q = 2 



(9) 
(10) 

(11) 



Although the sech' profiles can be adopted for any value of q 
between and 2 (e.g. de Grijs et al. 1997), it are these three 
models which are widely used in the literature to model the 
vertical structure in disc galaxies. For example, Schwarzkopf 
& Dettmar (2000) investigated the surface brightness of a 
set of 110 edge-on disc galaxies, and divided them into these 
three classes. They found that the intermediate sech-profile 
fitted the data (optical as well as near-infrared) best for 
55% of the galaxies, whereas the exponential and isothermal 
profiles fitted 35% and 10% of the galaxies. 

This would suggest to choose g = 1 as the typical 
value for our template model. However, in Section ^ we will 
demonstrate that the choice of q has only a minor eff'ect 
on the attenuation curve, such that we are relatively free 
to choose q to our convenience. We choose an exponential 
profile q — because it has some computational advantages 
(see Appendix 



3.3 The dust distribution 

Beside stars, we also need a dust component in the galaxy, 
which is determined by the opacity n{z). Again, the actual 
vertical distribution of dust in disc galaxies is difficult to in- 
vestigate. Direct information about the distribution of dust 
can be obtained by spatially resolved far-infrared or submm 
observations, where the dust emission is directly observed. 
During the last years, several nearby spiral galaxies have 
been modelled at these wavelengths using ISO and SCUBA 
data (Haas et al. 1998, Alton et al. 1998a, 1998b, Davies et 
al. 1999, Domingue et al. 1999, Trewhella et al. 2000). For 
NGC891, Alton et al. (2000) showed that the far-infrared 
emission in very good agreement with the dust distribution 
derived from radiative transfer calculations (Xilouris et al. 
1998), where an exponential vertical dust profile was as- 
sumed. 

It seems therefore obvious to assume that stars and dust 
have a similar distribution in the galaxy, but we allow the 
scaleheights to be different. We thus adopt a family of opac- 
ity functions 



k{z) 



qB{i'l)Ch'''^''{q(h) 



(12) 



which satifies the normalization condition (|^. For our tem- 
plate model we will thus adopt an exponential model for 
both stars and dust. The quantity ( represents the layering 
parameter, i.e. the ratio of the scaleheights of dust and stars 
(Disney et al. 1989). 

Combining the emissivity (^ and the opacity (|l^) we 
can calculate the source function S* (r) for our galaxy mod- 



To 



To 



(13) 



where St{tX,q) is normalized as 

/ s,{t,C,q)dt = l. (14) 
Jo 

The source function hence does not depend on either of the 
scaleheights of the stellar or dust distribution, but only on 
their ratio Therefore the same is true for the attenuation 
A{ii), which will thus depend on two geometry parameters, 
the shape parameter q and the layering parameter 

What is a representative value for the layering param- 
eter ? Normally, the interstellar matter sinks down to the 
central plane of a galaxy and forms an obscuring layer which 
is narrower than the stellar layer, such that < 1. This is 
observed in the Galaxy, where the scaleheights of the (thin) 
stellar disc and the dust disc are approximately 300 pc and 
100 pc respectively. In edge-on galaxies dust lanes suggest 
that also there the dust distribution is narrower than the 
stellar one. Using detailed radiative transfer modelling for 
seven such galaxies, Xilouris et al. (1999) find that the dust 
scaleheight is about half that of the stars. For our template 
model we set ( = 0.5, but in Section ^ we will consider a 
wide range in ^ in order to investigate the dependence of the 
relative star-dust geometry on the attenuation. 

The last parameter in our galaxy model is the total 
V band optical depth tv- Again, a realistic value for the 
optical depth in disc galaxies has been a subject of debate 
for a long time already. An extensive discussion can be found 
in e.g. Xilouris et al. (1997) and Kuchinski et al. (1998). 
The most widely supported idea is that spiral galaxies are 
transparant or moderately opaque if they are seen face-on, 
with a typical V band face-on optical depth of round unity. 
For our template model we adopt tv ~ 1, which corresponds 
to an optical depth of 3 in the far UV and 0.5 in the I band 
(Table |l|). 



4 THE INFLUENCE OF SCATTERING 
4.1 Approximations of the RTE 

The RTE as written in (^) is a partial integro-difl^erential 
equation, containing a differentiation along the path and an 
integration over the angle. The scattering term is the one 
that makes the RTE difficult, because this term is responsi- 
ble for the coupling of the RTE along different paths: due to 
the integration over p we cannot solve the RTE for different 
paths consecutively, but we have to solve it for all paths at 
the same time. Therefore it would be very practical if we 
could somehow get around the last term. 



4.1.1 No scattering 

One way is to neglect the scattered emission alltogether, 
i.e. to take into account the photons that are removed from 



* In Appendix ^ we explicitly calculate the source function for 
the exponential model. 
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the beam, but neglect the photons that are added to the 
beam due to scattering. This is the same as assuming that 
all the interactions between dust grains and photons are 
absorption events, and hence that there are no scatterings. 
We will denote the resulting attenuation curve as the ns 
attenuation. Mathematically this translates into setting uj — 
(the albedo is the ratio of the scattering efficiency to the 
total extinction efficiency), such that the last term in (^ 
vanishes, which turns the RTE into an ordinary differential 
equation 

/.|^(r,/i)=J(r,/i)-S.(r). (15) 

With the appropriate boundary conditions it can directly be 
solved, and the ns attenuation becomes 



A''\fi) = -2.5 log 



S,iT) exp ( -- ) dr. 



(16) 



A second way of getting around the scattering term is 
to neglect scattering completely, both the scattering absorp- 
tion and the scattering emission. The photons that would 
normally be scattered out of the beam, are thus assumed to 
continue on their path. For the RTE this means that not 
only the last term vanishes (the photons added to the beam 
due to scattering), but also the fraction of the first term (the 
photons removed from the beam due to scattering). In this 
case the RTE becomes 

M|^(r,^) = (l-c.)/(r,M)-5.(r). (17) 

This kind of interaction can also be described as completely 
forward scattering, because also then a scattering event has 
no net effect on a beam. The ARF for pure forward scatter- 
ing is ^'(/^, fJ.') ~ 2S(j£-- /i'), and inserting this in (^) we find 
the same equation (Oj). This equation is readily solved, and 
the fs attenuation is 



2.5 log / S^t) exp 



(1 



dr. 



(18) 



It is fairly straightforward to see that fs attenuation is in 
fact equivalent to ns scattering, but with the total optical 
depth replaced by an effective optical depth TcH = (1 ~uj)tq. 
If we thus obtain an expression for the ns attenuation, we 
immediately find the corresponding one for the fs attenua- 
tion. 



4- 1-2 (Approximate) isotropic scattering 

In radiative transfer problems, much attention has always 
been paid to isotropic scattering. In some physical situations 
the assumption of isotropy can be correct. For example, ac- 
cording to Bruzual et al. (1988) and Corradi et al. (1996), the 
asymmetry parameter is close to zero in the near-infrared K 
band, such that scattering is virtually isotropic therej^ An- 
other useful situation where the assumption of isotropy can 
be applied is a medium where photons have a large probabil- 
ity to be scattered more than once, because multiple scatter- 
ing tends to wash out the anisotropy effects. This principle 
has been applied by Vansevicius et al. (1997), who solve the 

t Note however that this assumption is not always followed, e.g. 
Gordon ot al. (1997) find gn = 0.43. 




Figure 3. The V band attenuation curves A{fi) as a function of 
H, corresponding to those in Figure ^. 



RTE in a disc galaxy using the iteration method. For the first 
scattering they adopt the general anisotropic scattering, for 
the next scattering events isotropic scattering is assumed. 

For isotropic scattering the distribution of the angles 
after a scattering event is uniform, such that the ARF is 
simply ^(/i,/i') = 1. A proper treatment of isotropic scat- 
tering hence does not change the integro-differential charac- 
ter of the RTE, and still requires the solution of the RTE for 
all path together. There are however two ways to approx- 
imate isotropic scattering such that the integration in the 
scattering term disappears. 

The first possible approximation to isotropic scattering 
is the so-called isotropic two-stream approximation, intro- 
duced by Code (1973)|. In the ts approximation a photon 
always remains on the same path, and after a scattering 
event it will either move in the same direction (i.e. forward 
scattering), or it moves in the opposite direction (i.e. reflec- 
tion). The probability of forward scattering and reflection 
are equal. The ARF corresponding to this kind of scattering 
is thus infinitely peaked in two directions, and can be writ- 
ten as vE'(/i, /i') = 5{fx — /i') -I- S{fi + fi'). The intensity for a 
certain value of jj, can then be found by solving the coupled 
set of RTFs 



dl. , 
dl 



(l-|) /(r,M)-|/(r,-/i)-5.(r) 



-M|^(r,-M)= (l-|)/(r,-/i)-|/(r,M)-S.(r) (19) 
For each inclination angle ji this set of ordinary differential 



t Code (1973) applied the two-stream approximation to calculate 
the flux escaping from a spherical dust envelope surrounding a 
star, and found a very good agreement with numerical results. He 
also proposes a second approximate solution, which works even 
better. However, this method cannot be applied to general plane- 
parallel (and other) geometries in a straightforward way, and is 
therefore not considered in this paper. 
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Figure 2. The various attenuation curves A{fi) as a function of wavelengtli. The left and right panel correspond to fi = 0.25 and fJ, = 1 
respectively. The adopted model is the template model. 



equations can be solved directly, such that we can calculate 
the two-stream attenuation curve A^''{fi). 

Another approximation for isotropic attenuation has 
been introduced by Natta & Panagia (1984), and after- 
wards adopted by various authors (e.g. Guiderdoni & Rocca- 
Volmerange 1987, Calzetti et al. 1994). Natta & Pana- 
gia (1984) suggest that a good approximation can be 
achieved by considering an attenuation equivalent to ns at- 
tenuation, but replacing the total optical depth by an ef- 
fective optical depth r^fr — \/l — uj to- The approximate 
isotropic attenuation or ai attenuation is thus 

A^'ifi) = -2.5 log S,(r) exp (^-^liH^j dr. (20) 

Although there is no physical principle behind this approx- 
imation, it is attractive due to its simplicity. 

4.2 Comparison of the attenuation curves 

In the previous section we have described five approxima- 
tions for the attenuation curve A{fi). These are compared 
in the Figures I and |, where we plot them respectively as 
a function of wavelength for two different angles /i, and as 
a function of angle for the V band. 



4-2.1 Scattering versus no scattering 

We first compare the attenuation curve with the ns approx- 
imation, i.e. the situation where the scattering emission is 
neglected. Qualitatively, the difference between the ns and 
all other attenuation curves is obvious: by not taking into 
account the scattered emission, the attenuation is overesti- 
mated in all directions. Because the probability of scattering 
is slightly higher than the probability of absorption at op- 
tical wavelengths (see Table |l|), it will be no surprise that 
this overestimation compared to the true attenuation can 



be considerable. In particular for face-on directions, we ob- 
tain a difference of nearly half a magnitude for our template 
model. Only for high inclinations, the ns approximation is 
more or less satisfactory. 

Another way to understand the effect of scattering is 
to compare the true attenuation curve with the fs attenu- 
ation curve, i.e. the situation where scattering is neglected 
completely. Particularly in Figure ^ the effect of including 
scattering on the attenuation curve is clearly observable: the 
attenuation becomes larger for small /i and smaller for large 
jj,. The overall effect of scattering is thus apparently that 
photons are removed from lines of sight with a high inclina- 
tion and sent into face-on directions, where they leave the 
galaxy. This is due to the fact that the optical depth along a 
path with inclination /i is proportional to l//i. Photons that 
initially (or after a scattering event) move on a path with 
large inclination have a large possibility to interact with a 
dust grain. If they are scattered into a direction which is 
nearly face-on, the probability to interact with another dust 
grain is much smaller, such that they can easily leave the 
galaxy. 

This overall effect of scatterig is responsible for a rather 
strange behaviour for face-on galaxies. If we look carefully at 
the right panel in Figure ^, we see that the true face-on at- 
tenuation becomes zero at the / band. The loss of radiation 
due to absorption (0.2 mag ~ 20 per cent) is thus com- 
pletely balanced by the gain due to scattering. This is not 
due to some particular properties of the dust at the I band, 
but only to the smaller optical depth. When we consider a 
smaller total optical depth, the extinction can even be neg- 
ative, i.e. a face-on galaxy can be brighter than it would be 
without dust. This is illustrated in Figure ^, where we plot 
the face-on attenuation in the l/-band as a function of the 
total optical depth. Using the iteration method such that 
the contribution of each partial intensity can be considered 
separately. The attenuation corresponding to the photons 
that leave the galaxy directly, i.e. the ns attenuation, is of 
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Figure 4. The V band face-on attenuation A{1) as a function 
of the total optical depth tq. Shown are the attenuation corre- 
sponding to photons that leave the galaxy directly, directly or 
after one scattering (direct -|- Is), directly or after one or two 
scatterings (direct -|- Is -|- 2s), and the total attenuation. Except 
for the varying total optical depth, the galaxy models have the 
same parameters as our template model (5 = and ^ = 0.5.) 



course positive everywhere. But adding the photons which 
leave the galaxy after one scattering event, the attenuation 
already is negative for tq < 0.25. The total attenuation re- 
mains negative for tq < 0.5, a result in correspondence with 
Bruzual et al. (1988) and Di Bartolomeo et al. (1995). 



4-2.2 Anisotropic versus isotropic scattering 

At first sight one would expect a significant difference be- 
tween the isotropic and the anisotropic attenuation curves. 
First, at optical and UV wavelengths the anisotropy pa- 
rameter is fairly large {g ~ 0.5), and for these wavelengths 
the Henyey-Greenstein phase function is far from isotropic 
(e.g. Bianchi et al. 1996, Figure 2). Second, the total opti- 
cal depth in our galaxy model is only moderate, such that 
the anisotropy effects should not be washed out by multiple 
scattering. However, Figure ^ and ^ show that the isotropic 
attenuation approximates the anisotropic attenuation very 
well, with differences between both attenuation curves only 
a few hundredths of a magnitude. 

In Figure Al of Paper I we plotted the Henyey- 
Greenstein ARF for various values of /i'. For nearly face-on 
lines-of-sight n' , the distribution of the new angles /j, is far 
from uniform. For very inclined lines-of-sight however, the 
ARF is fairly uniform, due to the averaging out over all az- 
imuths. For photons which are scattered from such lines-of- 
sight, the scattering seems thus approximately isotropic. In 
our galaxy model, the majority of the scattering events oc- 
curs on paths which are nearly edge-on, because along such 
paths the optical depth is much larger than on paths which 
are nearly face-on. This explains why even for the moder- 
ate optical depth of our galaxy, the differences between the 
isotropic and anisotropic attenuation are fairly small. 




- 1.5 



- 1 



- 0.5 



Figure 5. The source function 5, (t) of the family of sech' mod- 
els as a function of the parameters q and (,'. In each of the panels 
the template model is the black curves, and we vary one of the 
geometry parameters. Left: variation of the shape parameter q, 
taking the values 0, 1 and 2. Right: variation of the layering pa- 
rameter taking the values 0.25, 0.5, 1 and 1.5. 



As mentioned, there is little gain in the observation 
that the isotropic attenuation curve approximates the true 
anisotropic one so well, because the approximation does not 
simplify the RTE. However, if the ts or ai attenuation curves 
were to be satisfying approximations for the isotropic at- 
tenuation curve, the operational gain would be significant. 
Unfortunately, the Figures |^ and g show that the approxi- 
mations are not satisfying. The two-stream approximation 
is both qualitatively and quantitatively comparable to the 
fs approximation, i.e. an overestimation of the attenuation 
in face-on directions of more than 0.1 mag, and an un- 
derestimation of the same order for highly inclined direc- 
tions. The ai attenuation curve lies between the ns and fs 
attenuation curve, which is logical if we compare the ex- 
pressions ([1^, dlq ) and ([20[). This means that the approx- 
imation is satisfying for highly inclined directions, but the 
attenuation is severily overestimated in the face-on direction 
(A^« 0.25 mag). 



5 THE INFLUENCE OF GEOMETRY 

Our models contain two parameters which specifically deter- 
mine the geometry of the galaxy model: the shape parameter 
q determines the actual distribution of stars and dust near 
the galaxy plane, whereas the layering parameter determines 
the relative distribution of stars and dust. In this section we 
investigate to which degree these geometry parameters de- 
termine the attenuation curve of the galaxy. 

5.1 The shape parameter q 

The RTE (^ shows that the source function completely de- 
termines the attenuation curve, i.e. once S,{t) is known, 
we can calculate A(/i). Therefore, if we want to know the 
influence of q on the attenuation curve it is interesting to 
firstly study its influence on the source function. In Figure ^ 
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we plot S,{t) for various models with different shape pa- 
rameters. The source function depends only weakly on the 
parameter q, more precisely the shape of 5** becomes slightly 
rounder at the centre of the galaxy as q increases. It seems 
logical to predict a similar insensitivity on q for the atten- 
uation curves. This is confirmed by our calculations. Even 
the differences between the attenuation curves of the two ex- 
treme cases, the exponential disc and the isothermal sheet, 
are almost negligibly small. 



5.2 The layering parameter 

The layering parameter is the ratio of the scaleheights of 
dust and stars and hence determines the relative distribution 
of the two components. Variation of this parameter yields a 
wide range of dust-stars geometries. This can be seen in 
Figure |5|, where the source function is plotted for various 
models with different layering parameter i^. According to 
the shape of the source function for a particular value of ^, 
we can classify the different geometries into three classes. 

For C = li the stars and dust have the same scaleheight 
and thus the same spatial distribution. The source function 
is then independent of the depth in the galaxy, and due 
to the normalization (^) we find S{t) = I/tq. In general, 
each galaxy model with a similar distribution of stars and 
dust will have such a constant source function. Because the 
attenuation curve is completely determined by the source 
function only, all galaxy models with a similar distribution 
of stars and dust will be equivalent. In particular, our model 
will be equivalent to the homogeneous slab, i.e. a finite galaxy 
model with stars and dust homogeneously mixed. Such a 
model has been considered by many authors (e.g. Bruzual 
et al. 1988, Calzetti et al. 1994, Di Bartolomeo et al. 1995) 
as an approximation of a disc galaxy. 

For C < 1 the dust has a smaller scaleheight than the 
stars, such that the attenuation will predominantly occur in 
the central regions. As mentioned in Section 3.c , this is the 
most realistic range for the layering parameter. If < 1, 



the source function has its minimal value in the centre of 
the galaxy, and diverges at the edges. If ( tends to zero, 
the dust will form an infinitely thin obscuring layer in the 
central plane of the galaxy. We will refer to this model as 
the thin dust disc model. 

For C > 1 the dust has a larger scaleheight than the 
stars, and the attenuation will thus occur relatively more 
in the outer regions of the galaxy. The source function is 
then maximal in the centre of the galaxy, and vanishes at 
the edges of the galaxy. If grows larger, the bulk of the 
dust attenuation will occur further away from the central 
plane of the galaxy. In the limit — > oo the dust is infinitely 
thinly distributed over all space, and will effectively behave 
as a set of two obscuring layers on both sides of the galaxy. 
This geometry, where the obscuring material is located only 
between the observer and the source, is known as the over- 
lying screen model. It is the analogue of attenuation of stars 
in the Galaxy. Various authors advise against the use of this 
model (Disney et al. 1989, Witt et al. 1992). 

The source function is thus strongly dependent on the 
relative distribution of stars and dust, whereas the actual 
shape of the distribution is only of minor importance. We 
will now investigate in detail is whether this also means that 



the influence of the layering parameter on the attenuation 
curve is significant. 

5.2.1 Effects on the fs attenuation curve 

First we will only consider true absorption and neglect scat- 
tering, and we will hence investigate the influence of C, on 
the fs attenuation curve. For the exponential model q = 0, 
the source function is fairly simple, and the integration ( |l8| ) 
can be performed analytically. One flnds 

Af=(^) . 0.543 (i^^ - 2.5 log Wc 



(21) 



This expression is derived in Appendix ^ and the definition 
and some properties of the function (x) are described in 
Appendix ^ The advantage of the expression ( pl| ) is that it 
easily allows us to study the attenuation curve as a function 
of (. For every value of x, W^(a;) is a decreasing function of 
moderately decreasing for small x and strongly decreas- 
ing for large x. This implies that j4^°(/x) will be a increasing 
function of (,. The more extended the dust distribution, rel- 
ative to the stellar one, the larger the attenuation and thus 
the fainter the galaxy appears. In particular, the overlying 
screen model is the most efficient and the thin dust disc the 
least efficient absorbing geometries. This behaviour is illus- 
trated in Figure where we plot the fs attenuation curve 
A^^(^j,) as a function of fi for a set of layering parameters 
For highly-inclined directions (small values of /i) the fs at- 
tenuation changes considerably as a function of the layering 
parameter. 

5.2.2 Effects on the true attenuation curve 

Let us now investigate the influence of on the true at- 
tenuation curve, which has to be calculated numerically. In 
Figure ^ the results are shown: we see a different behaviour 
between the fs and the true attenuation curves. Whereas the 
fs attenuation is a increasing function of ( in all directions, 
the true attenuation is not: for small inclination angles, the 
true attenuation decreases as a function of (. This difference 
is also illustrated in Figure |^, where we plot the fs and true 
attenuation explicitly as a function of the layering param- 
eter for different inclination angles. The fs curves all rise, 
whereas the true attenuation curve rises for small values of 
H, but decreases in the face-on direction jj, — 1. 

If scattering is not taken into account, galaxies with 
an extended dust distribution are more efficiently obscured 
than galaxies with their dust concentrated in the central 
region. However, if scattering is taken into account, the ef- 
ficiency turns around for face-on directions. In face-on di- 
rections, galaxies with their dust concentrated in the cen- 
tral region are more efficiently obscured than galaxies with 
an extended dust distribution. How can this be understood 
physically? In Section ^ we showed that scattering biased 
the preferential direction for photons to leave the galaxy: 
it is easier to leave it face-on than for highly-inclined di- 
rections. In particular photons which are scattered near the 
edge of the galaxy are responsible for this: if a photon is 
scattered there into a face-on direction, its probability to 
leave the galaxy without another interaction is high. Now, 
the more extended the dust distribution in the galaxy, the 
more dust in the outer regions, and the more photons which 
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Figure 6. (a.) The fs attenuation AY(fi) curve as a function of /i for some different values of the layering parameter which are 
indicated. The black curves denoted by = and (J = oo correspond to the thin dust disc and overlying sheet models respectively, 
(b.) Same as a., but now the true attenuation curve, (c.) The true and fs attenuation curves as a function of the layering 

parameter C^. The solid lines represent the true attenuation curve, the dotted lines are the fs attenuation curves. The results are shown 
for four angles, which are indicated. 



are scattered out of the galaxy in nearly face-on directions. 
Relative to galaxies with concentrated dust, galaxies with 
an extended dust distribution on the one hand efficiently 
absorbe photons and on the other hand they scatter a con- 
siderable amount of photons in the face-on direction. 

Which of these two effect dominates is determined by 
the total optical depth. For small optical depths the scat- 
tering effect will be more important, and A{fi) will decrease 
with increasing for nearly-face-on directions. In contrary, 
for large optical depths, the absorption will dominate and 
will be an increasing function of ^. We can thus con- 
clude that 

• inclined galaxies with a moderate optical depth are 
most efficiently obscured if the dust is more extended than 
the stars. 

• nearly face-on galaxies with a moderate optical depth 
are most efficiently obscured if the dust is concentrated in 
the central regions. 

• optically thick galaxies are most efficiently obscured if 
the dust is more extended than the stars. 

Notice that the difference between small and large optical 
depth can also be considered for one single galaxy, when 
we study the attenuation as a function of wavelength. This 
is shown in Figure where we plot the attenuation curve 
Eis a function of wavelength for two angles p and for three 
models with a different layering parameter. For a total V 
band optical depth tv ~ 1, the attenuation behaves as an 
optically thick galaxy at the blue side of the spectrum, and 
as an optically thin galaxy in the red wavelength range, with 
the crossover between the two regimes occurs around 300 nm 
(ro ^ 2). 

These results can be compared with those found by 
Bruzual et al. (1988), who adopted a sandwich model (i.e. a 
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Figure 7. Three attenuation curves A(fj,) as a function of wave- 
length for three models with a different layering parameter. The 
values for f are 0.1, 0.5 and 1, the other parameters are those of 
the template model {ry = 1 and q = 0). The upper panel and 
the lower panel correspond to ^ = 0.25 and fi = I respectively. 
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slab where stars and dust are homogeneously mixed, sur- 
rounded on both sides by a slab containing stars only) . They 
varied the relative thicknesses of the slabs, and found a sim- 
ilar result of the dependence of the attenuation curve on 
the layering parameter (only shown for the face-on direc- 
tion /i = 1). Bruzual et al. (1988) concluded that the differ- 
ences in the attenuation curve between the various models 
have the same order of magnitude than the uncertainties 
due to an imprecise knowledge of the optical parameters of 
the dust, and hence that a homogeneous slab is a satisfy- 
ing model for disc galaxies. We agree that in the face-on 
direction, quantitatively the difference in attenuation be- 
tween centrally concentrated and extended distributions is 
of the order A A ~ 0.1 mag for optical wavelengths. These 
differences become even smaller if we constrain the range of 
^ to realistic values. The layering parameters obtained by 
Xilouris et al. (1999) all lie between 0.30 and 0.75 at optical 
wavelengths, and for these values AA{n) is only of the order 
of a few hundredths mag. Also for moderately inclined galax- 
ies, the attenuation curve will not seriously depend on ^ if 
it this parameter is constrained to a realistic range. Only if 
the optical depth is considerably higher, due to (a combina- 
tion of) a larger amount of dust, a more opaque wavelength 
region or a large inclination, the uncertainty of the relative 
geometry of dust and stars causes a considerable uncertainty 
on the attenuation. 



6 DISCUSSION AND CONCLUSION 

The aim of this paper was to investigate how the RTE can be 
simplified and what errors are introduced by doing so. More 
than a decade ago already, Bruzual et al. (1988) investigated 
to the effects of scattering on the attenuation in disc galaxies 
(albeit for a simple plane-parallel model), and Disney et 
al. (1989) showed that the adopted geometry of the galaxy 
model strongly influences the attenuation (without taking 
scattering into account however). We wanted to reexamine 
these effects in detail and in the same manner, i.e. within 
a realistic set of disc galaxy models and properly taking 
multiple scattering into account. 

In Section ^ we showed that the nature of the physi- 
cal process of scattering is such that it connects different 
lines of sight with each other. The overall effect of scat- 
tering is that many photons which originally move on a in- 
clined line-of-sight, leave the galaxy in the face-on direction. 
Any approximate method to simplify the RTE where pho- 
tons cannot leave their initial path (such as the ns, fs, ts 
and ai approximations described in Section will not be 
able to reproduce this behaviour. In particular for face-on 
galaxies, the optical depth will always be underestimated 
by such approximations, as illustrated in Figure |^. If we 
compare the attenuation curves corresponding to isotropic 
and anisotropic scattering however, we see that the differ- 
ence between them are only of the order of a hudredth of 
a magnitude. Anisotropy effects are essentially washed out. 
Inversely it will be impossible to determine an ARF from 
photometry. 

The effects of scattering are also important if we in- 
vestigate the influence of the star-dust geometry on the at- 
tenuation in disc galaxies. If scattering is not taken into 
account, as in Disney et al. (1989), we find that the atten- 



uation is more effective if the dust distribution is extended. 
For a same amount of dust, the larger the layering param- 
eter, the larger the attenuation. In particular this means 
that the overlying screen geometry is the most efficient ob- 
scuring geometry, whereas galaxies with a dust scaleheight 
smaller than the stellar scaleheight, are less efficiently ob- 
scured. However, when scattering is taken into account, this 
picture does not hold anymore. We find the remarkable ef- 
fect of scattering that in the face-on direction, galaxies with 
a moderate optical depth (tq < 2) are most efficiently ob- 
scured by a centrally concentrated dust distribution. This is 
very interesting because the optical depth of spiral galaxies 
is thought to be of the order unity in optical bands. 

The effects of scattering and geometry not only need to 
be investigated in a qualitative way, it is important to know 
them quantitatively. This allows us to answer questions as 
how large our error is when we approximate anisotropic scat- 
tering by assuming it is forward, or when we assume a dust 
scaleheight equal to that of the stars, where in reality it is 
only half of it. 

The magnitude of these errors depends on the incli- 
nation. In the face-on direction, the errors induced by not 
properly taking scattering into account are largest. In the 
ns approximation, the induced errors rise up to nearly half 
a magnitude, for forward scattering they are of the order 
of 0.15 mag at optical wavelengths (Figure ^). The errors 
introduced by uncertainty about the star-dust geometry are 
only of the order of a few hundredths of a magnitude (Fig- 
ure ^), and thus significantly smaller. For inclined direc- 
tions, the optical depth along lines-of-sight towards the ob- 
server is much larger. The attenuation curve then becomes 
less dependent on the angular redistribution function, in 
agreement with Di Bartolomeo et al. (1995). The differences 
between the attenuation curves corresponding to different 
star-dust geometries grows significantly. For example, for a 
galaxy with an inclination of 75° (/i — 0.25), the differ- 
ences between the true, isotropic, forward and two-stream 
attenuation curves become about 0.1 mag at optical wave- 
lengths, and only a few hundredth of a magnitude in the UV 
(Figure ^a). However, assuming a homogeneous slab model 
instead of our template model, the induced error on the at- 
tenuation curve is about 0.3 magnitudes at optical and more 
than half a magnitude at UV wavelengths. For face-on galax- 
ies it is thus important to include scattering in a proper way 
in the modelling, for inclined galaxies it is important to have 
a fairly good idea of the star-dust geometry of the galaxy. 

In the light of the adopted models, two remarks are ap- 
propriate. (1) By working in plane-parallel geometry we have 
the advantage that we could use four independent methods 
to solve the RTE, such that the accuracy of the obtained 
results can be checked carefully. We are well aware of the 
fact that plane-parallel galaxy models do not represent very 
realistic galactic discs. More detailed modelling requires an 
exponential fall-off in the radial direction (Freeman 1970), 
a distinction in attenuation between arm and interarm re- 
gions (e.g. White et al. 2000), and an inclusion of the effects 
of clumping of both stars and dust (Witt & Gordon 1996, 
2000; Bianchi et al. 2000). Nevertheless we are convinced 
that our modelling can contribute to the understanding of 
the mechanism of radiative transfer, and we have all reasons 
to believe that our results would hold qualitatively in more 
complex models. (2) We know that scattering and geometry 
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are not the only two elements that introduce uncertainties 
to the attenuation of disc gala:xies. In Article I we solved 
the RTE for identical galaxy models, but with two different 
sets of optical dust parameters found in the literature. We 
obtained differences between the attenuation curves of the 
order of a few tenths of a magnitude. Uncertainty about the 
optical properties of the dust hence also introduces errors 
to the attenuation curve. However, these errors are inde- 
pendent of those introduced by an inaccurate treatment of 
scattering or geometry, and they hence do not affect the 
conclusion of this work. 

As a final remark, we re-iterate a warning that has also 
been issued by others: do not use too simplistic models to 
determine the opacity in disc galaxies. There are many ways 
to investigate the amount of attenuation or the optical thick- 
ness of disc galaxies. Following Xilouris et al. (1997), we can 
divide them into two categories, the small N approach and 
the large N approach. The small A'" approach tries to de- 
termine the opacity of individual nearby disc galaxies by 
constructing detailed radiative transfer models. Obviously, 
these methods would reveal wrong results if scattering is ne- 
glected. The large N approach aims at determining a mean 
value for the optical depth in disc galaxies, by statistically 
studying a large sample of galaxies. Of such a sample of 
galaxies a certain observable is then compared to a galaxy 
model, in order to derive the optical depth. However, these 
methods are difficult to interpret because of selection effects 
(Davies et al. 1993), and there arc often model-dependent. 
For example, Saunders et al. (1990) estimate an average 
effective optical depth of disc galaxies by fitting a screen 
model to the optical and far-infrared luminosity functions. 
Trewhella et al. (1997) redid these calculations with a sand- 
wich model, and found a very different conclusion. We want 
to stress that also in such studies the effects of scattering 
need to be taken into account, because they are, except for 
very inclined galaxies, at least as important as the effects 
of geometry. An often used prescript to take scattering into 
account, is to determine an effective optical depth T^ff for 
a galaxy model without scattering and compensate the ob- 
tained value for scattering. These compensations are simply 
multiplications, such as tq — Tcft/{l—uj) or to = Tcb/^/I — ui, 
which can but not necessarily need to have a physical justifi- 
cation. With our paper we hope to have clearly showed that 
such approximations are not able to accurately reproduce 
the effects of scattering. We therefore argue that scatter- 
ing needs to be accounted for properly in statistical studies 
concerning the opacity of disc galaxies. 
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APPENDIX A: THE EXPONENTIAL MODEL 

In this appendix we take a closer look at the exponential 
model q = 0, characterized by 
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Disney et al. (1989) already considered a similar galaxy 
model (their Triplex model). They managed to calculate 
an expression for the observed intensity if scattering is not 
taken into account. We will derive the solution of the RTE 
by direct integration of the formal solution (p^. Therefore 
we first need an explicit expression for the source func tion 
S*(t). The optical depth is obtained by combining ( A2) 
with (|), 



To 

2 

To 



and inverting this relation we obtain 



for 2 
for jz > 0, 



(A3) 



z{t) 



^(,h\n[2{l - t)] 



for t ^ i 
for i ^ t ^ 1, 



(A4) 



where t = t/tq is t he rela tive optical depth. Introducing this 
expression in (Al) and (A2) we find an expression for the 



stellar source function S*(t), 



^(2t)^-i 

To 

i-[2(l-t)] 



for s; t ^ ^ 



for i s; t ^ 1. 



(A5) 



The source function is of the form (|13|), and the normaliza- 
tion condition ( p^ ) is satisfied. The fact that 5'*(t) has a 
fairly simple power law form, will allow us to calculate an 
expression for the ns attenuation, by direct integration of 
the formal solution ([l6|. Inserting S, (t) in ( |l^ ) we obtain 



A''%n) = 0.543 



To 



2.5 log C / X 

10 



cosh 



(x-1) 



Ax, (A6) 



which becomes, using the function VV{(a;) defined as in Ap- 
pendix B, 



yl"=(M) = 0.543 — -^2.5 log Wcf — 



(A7) 



This result is in agreement with the solution obtained by 
Disney et al. (1989). Knowing this expression, we immedi- 
ately also have an expression for the fs and ai attenuation 
curves 



A'^itJi) = 0.543 ^"^ - 2.5 log Wc 



(1 — uj)to 



2fi 



(A8) 



A-i(^) = 0.543 _ 2.5 log Wc ( yll^) . (A9) 
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Figure Bl. The function yV/^ix) as a function of x for different 
values of fj. 



APPENDIX B: THE FUNCTION Wc(a;) 

For C > and for all values of x, we define the function 
Wc(a;) by the integral 



>Vc(a;)=C/ (1-t)^" cosh(2:f)dt 



(Bl) 



For natural values of the integral can immediately be 
solved in terms of elementary functions, for example 

Wi(x) = ^ (B2) 



W2(x) 

W3(x) = 6 



X 

cosh X 



sinh X -\- X 



(B3) 
(B4) 



For non-integer values of C, the integral can be calculated by 
expanding the cosh function. 



14- 



r(C + i) 
r(c + 2fc + i) 



(C + i)(C + 2) 



+ 



(B5) 



(B6) 



(C+l)(C + 2)(C + 3)(C + 4) 

Using d'Alembert's convergence theorem one can easily 
show that this power series converges absolutely. It can be 
written as a hypergeometric function, 

„2 ^ 



Wc(a;) = 1F2 1 



C + 1 C + 2 X' 



(B7) 



Considering the expansion (B5) we can also consider the 
special cases C = and ^ = 00, 



Wo(a;) = cosh a; 
Woo (a;) = 1. 



(B8) 
(B9) 
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In Figure Bl we plot W((x) as a function of a; for a few values 
of It is an even function of x, increasing for positive x, 
the minimal value being W^(0) = 1 for all Regarded as a 
function of ( for fixed values of x, Wij(a;) it is a decreasing 
function, depending weakly on for small x and strongly for 
large x. 

This paper has been produced using the Royal Astronomical 
Society /Blackwell Science I^TJ^X style file. 
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